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J> , Abstract 

We calculate m„d = (m„ + md)/2, ms, and Jk in the quenched 
continuum limit with UV- filtered overlap fermions. We see rather small 
scaling violations on lattices as coarse as ~ 1 GeV and conjecture 
that similar advantages would be manifest in unqucnched studies. 



^ ■ 1 Introduction 

oo 

' Overlap fermions ^ satisfy the Ginsparg- Wilson relation (21 

00 ; 

ToD + Di5 = 0, = 75(1 -il)) (1) 

O ; P 

• which implies exact chiral symmetry at finite lattice spacing Pl ^^ . While 

theoretically clean, calculations with overlap fermions are considered a com- 
Q^l putational challenge. In a recent investigation it has been conjectured 

that a UV-filtered (smeared) Wilson kernel operator, as previously suggested 
in [01 [71 IHl ini 5 might substantially reduce the computational cost associated 
with obtaining continuum physics [10^. In ^ the focus was on technical 
^ ' aspects, but it is clear that the point relevant in practical applications is 

. whether UV-filtered ( "thick link" ) overlap fermions would extend the scaling 

region to significantly coarser lattices, even if one refrains from (/3-specific) 
tuning and sets p=l- The present note addresses this question by study- 
ing the continuum limit of the pseudoscalar masses and decay constants 
with quenched UV-filtered overlap quarks. No systematic comparison to 
plain ("thin link") overlap fermions is made, because we could not obtain 
reasonable signals for plain overlap fermions with p = 1 on our coarser lat- 
tices. Given the exploratory nature of the present investigation, we restrict 
ourselves to pseudoscalar meson correlators in the p-regime of Chiral Pertur- 
bation Theory (XPT) Our main results are the values of the quenched 
strange quark mass and K decay constant in the continuum 

m,(MS, 2 GeV) = 119(10)(7) MeV , fx = 170(10)(2) MeV (2) 
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as well as the quark mass and decay constant ratios 



rriud 



23.3(7.1)(4.5) , fK/fn = 1.17(4)(2) (3) 



where mud = (nT'u+nT'd) and no numerical input from XPT has been used. 
Throughout this note the first error is statistical and the second systematic, 
but the latter does not include any estimate of the quenching effect. 

2 Technical setup 

We use the Wilson gauge action. The massless overlap operator is P 

Dov = + - p){{dI - p){D^ - p))-^'''] (4) 

with the massless Wilson operator. The UV-filtered overlap is construc- 
ted by evaluating the Wilson operator on APE ^2 or HYP |13Lll4j smeared 
gauge configurations, resulting in an 0{a?) redefinition of the fermion action 
[T^ Ej. We use 1 or 3 iterations with smearing parameters oape = 0.5 or 
«HYP = (0.75, 0.6, 0.3) and the shift parameter p = l\s kept fixed. Based on 
(jlj the massive operator is defined through 

-Dov.m = V'^ ~ "2/7/ ■ ^' 

We set the lattice spacing with the Sommer parameter ^S], i.e. we give 
all intermediate results in appropriate powers of tq to facilitate comparison 
with other quenched studies. Only the final result will be converted into 
physical units assuming a standard value for tq. We have data at 4 different 
lattice spacings in matched boxes of physical volume x 2L with L~3ro. 
The couplings were chosen with the interpolation formula given in jl7j . The 
details of the simulation are summarized in Tab.^ 

For each coupling and filtering level we use 4 bare quark masses. Ideally 
one would choose them such as to always obtain the same 4 renormalized 
quark masses in tq units (or the same 4 pseudoscalar masses), but for this 
one would have to know the renormalization factor Zm = -^5 beforehand. 
Our values as summarized in Tab. [2 are not bad; our renormalized masses 
are roughly in the region \mP^^^ . . .m^^'^^. 





8^x16 


10^ X 20 


12^x24 


16^x32 


/? 


5.66 


5.76 


5.84 


6.0 


a[fm] 


0.188 


0.149 


0.125 


0.093 


#conf. 


30 


30 


30 


30 



Table 1: Simulation parameters and statistics. The box volume in physical 
units xT= (1.5 fm)'^ X 3.0fm (based on ro = 0.5fm) is kept fixed. 
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8^ 


X 16 


10^ 


X 20 


12^ 


X 24 


16^ 


X 32 


none 


0.16 


... 0.4 


0.16 . 


.. 0.4 


0.16 . 


.. 0.4 


0.08 . 


.. 0.2 


1 APE 


0.08 


... 0.2 


0.08 . 


.. 0.2 


0.04 . 


.. 0.1 


0.04 . 


.. 0.1 


3 APE 


0.04 


... 0.1 


0.04 . 


.. 0.1 


0.03 . 


.. 0.075 


0.03 . 


.. 0.075 


IHYP 


0.08 


... 0.2 


0.08 . 


.. 0.2 


0.04 . 


.. 0.1 


0.04 . 


.. 0.1 


3 HYP 


0.04 


... 0.1 


0.04 . 


.. 0.1 


0.03 . 


.. 0.075 


0.03 . 


.. 0.075 



Table 2: The 4 regularly spaced bare quark masses per coupling and filtering. 

In the course of this calculation we will need both Zs = Zp (to determine 
the renormalized quark masses) and Zy = Za (for the decay constants), 
where the alleged identity is specific for the massless overlap operator. 

To compute the scalar and pseudoscalar renormalization constant we 
follow the non-perturbative RI-MOM procedure as defined in ^5] and first 
applied to overlap fermions in jElEOl- We compute both Zs{n) and Zp{ii); 
in the latter l/m term is used to extrapolate the result to the chiral 

limit. It turns out that the values obtained are in good agreement even on 
our coarsest lattice. We fit the scalar renormalization constant to the form 

^Ri-MOM(^) = U{^,)Zf''' + const (a/x)2 (6) 

where U (fi) is the 4-loop running in the RI-MOM scheme as given in |2J 
and the second term is introduced to account for discretization effects. In 
other words the scalar renormalization constant after dividing out the 4-loop 
perturbative running should be flat, up to discretization effects, and Fig.^ 
shows that the latter are indeed non-negligible. The wiggles in the data 
signal rotational symmetry breaking on the lattice. We checked that within 
errors the slope disappears in proportion to (a/ro)'^. The phenomenological 
analysis below is based on ^^^"^^^^^(2 GeV), where /i = 2 GeV is realized 
through firQ = 5.06773. The systematic error is estimated by varying the 




Figure 1: Z^^-^^^^{n)/U{fi) on our coarsest and finest {(3 = 5.66, 6.0) lattice 
using the 1 HYP overlap operator with p=l. The solid line indicates a linear 
fit with range 3</iro<5 and 5</iro<9, respectively. 
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(3 = 5.66 


/? = 5.76 


13 = 5.84 


(3 = 6.0 


none 


ill-def. 


ill-def. 


6.14(15)(38) 


2.54(6)(12) 


lAPE 


3.05(7)(34) 


1.79(4)(6) 


1.60(3)(8) 


1.26(3)(8) 


3 APE 


2.05(7)(43) 


1.35(4)(8) 


1.25(2)(8) 


1.04(2)(7) 


IHYP 


1.71(4)(24) 


1.24(3)(6) 


1.23(2)(6) 


1.03(2)(7) 


3 HYP 


1.57(6)(25) 


1.16(3)(5) 


1.10(2)(7) 


0.98(2)(6) 



Table 3: Renormalization constant {2GeV) for Wilson glue and /0 = 1. 
The main effect of filtering is to bring it much closer to its tree- level value 1. 



fit range, by including additional l/p^ terms into the fit, and by comparing 
to the Zp^~^^^(2 GeV) data. A summary of our results, after conversion 
to (MS,2GeV) conventions, is presented in Tab.Ol Choosing a fixed p > 1 
would delay the breakdown of the unfiltered version - see [H] for details. Note 
that the Zs factors of all UV-filtered operators are much closer to 1, even 
when compared to the unfiltered overlap operator with tuned p |19[ I22| I23j . 
This suggests that one should be able to compute renormalization constants 
perturbatively, as was done in [2] in a slightly different setup. 

The second ingredient is the axial- vector renormalization constant Za- 
Here we use the values given in (coming from a PCAC renormalization 
condition), complemented by a /9 = 5.76 column included in |24j. It turns out 
that the values in this column are in fair agreement with the prediction by 
the Pade curve given in ;5,, which is another indication that for the filtered 
overlap operator lattice perturbation theory might work rather well. 



3 Physical results 

To extract meson masses and decay constants we compute the correlators 

CT„T,{t) = ^(v;i(o)ri^2(o)^2(x,t)r2V'i(x,t)) (7) 

X 

where 

; 1 + 7575 , 

^ = ^ V' (8) 

denotes the "chirally rotated" quark field Hj. Specifically, we consider 

(S5)(t)=Ci,i(t), (PP)(t) = C,,,^3(t), (ylo^o)(t) = C^57o,757oW (9) 

and we extract the mass and decay constant of the pseudoscalar meson from 
fits to the (^o^o) and {PP)-{SS) channels. We do not consider (PP), 
since it is known to be contaminated by zero mode contributions |2()| I26j . 

Generically, we see plateaus in the effective mass from about Nt/4: on, 
with the exception of the unfiltered operator, where the plateau sets in later 
and is less pronounced. Our central values stem from a fit to the (AqAq) 
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Figure 2: Mprg versus the bare quark mass (777,1+771.2)^0 on our coarsest and 
finest (/3 = 5.66, 6.0) lattice for the IHYP operator. The masses come from 
a fit to the (AqAq) correlator in the range 5A'7^/16<t<llA'2^/16. The solid 
curve represents a fit to the functional form 1)10(1 and the horizontal lines 
indicate the physical M.^rQ and M^tq values. The latter are used to read 
off the fitted 2mudfQ and {ms+mud)ro-, respectively. 

correlator in the range 5A^t'/16 < t < llA^r/16. The theoretical error is 
dominated by the comparison to {PP) — {SS). Including only variations of 
the fit range (and below the error on Zs), it would be much smaller. 

In Fig. [21 the pseudoscalar meson mass squared is plotted versus the bare 
quark mass. Using regularly spaced quark masses, each 7711+7772 combination 
is realized in several ways and the pertinent Mp are in excellent agreement. 
In other words, isospin breaking effects are completely negligible. Therefore 
we fit the quark mass dependence of the pseudoscalar meson mass with the 
resummed quenched XPT expression [^Tj 

M|, = A (7771 +7772) ^/(^+'^) + B (7771+7772)^ (10) 

which is strictly true only for degenerate masses. Our 6 fiuctuates wildly 
and is thus trivially consistent with 0.2 

The next step requires some experimental input. Since we do not see 
any isospin breaking effects and electromagnetic corrections are of the order 
of 1%, while our statistical errors turn out to be roughly 10%, we decided 
to use the charged meson masses as input. Sticking to the identification 
ro = 0.5fm we use M/^ro = 493.7 MeV- 0.5 fm = 1.251 and M^ro = 139.6 MeV- 
0.5 fm = 0.3537 as input. This gives the bare (7771 +7772)7'o values indicated 
by a horizontal error bar in Fig. [3 Multiplying the bare masses with the 
Zm = Z^^ obtained earlier, we extract the renormalized quark masses. Our 
results for nis+niud and 2777^^ are reported in Tabs. 0] and respectively. 

To finally push to the continuum we perform a combined fit (including 
all couplings and smearing levels) with a common continuum value and 
individual 0{a?) terms. We do this for the quantities {rns+mud)fo-, 2mudfo 





f3 = 5.66 


f3 = 5.76 


P = 5.84 


13 = 6.0 


none 
1 APE 
3 APE 
IHYP 
3 HYP 


ill-def. 
0.244(34) (45) 
0.260(21)(47) 
0.301(25)(39) 
0.281(20)(40) 


ill-def. 
0.366(40)(50) 
0.311(32)(43) 
0.336(35)(67) 
0.321(31)(24) 


0.219(34)(181) 
0.260(33)(19) 
0.223(32)(20) 
0.240(28)(23) 
0.247(31)(21) 


0.378(36)(28) 
0.329(27)(31) 
0.312(26)(26) 
0.321(26)(28) 
0.321(26)(22) 



Table 4: (m<j + m„(i)ro in (MS,2GeV) conventions. For the p = l thin link 
action at /? = 5.66, 5.76 no Zjn = Zg is available (cf. Tab. 3). 





(3 = 5.66 


l3 = 5.76 


(3 = 5.84 


(3 = 6.0 


none 
lAPE 
3 APE 
IHYP 
3 HYP 


ill-def. 
0.006(07)(09) 
0.012(05)(17) 
0.021(07)(05) 
0.014(08)(05) 


ill-def. 
0.017(12)(13) 
0.014(10)(18) 
0.021(08)(09) 
0.028(12)(09) 


0.006(05)(08) 
0.007(08)(04) 
0.007(07)(01) 
0.011(07)(02) 
0.009(06)(05) 


0.027(09)(03) 
0.027(07)(04) 
0.023(06)(03) 
0.026(07) (04) 
0.022(05) (03) 



Table 5: 2mudfQ in (MS,2GeV) conventions. For the p = l thin link action 
at /? = 5.66, 5.76 no Zm = Zg is available (cf. Tab. 3). 




Figure 3: Renormalized (MS, 2 GeV) masses (ms+m„d)ro and 2mudfo versus 
{a/r^f'. Errors are statistical, the 4 filterings are correlated (same configs.). 
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Figure 4: Pseudoscalar decay constant /pro versus (Mpro)^ on our coars- 
est and finest (/? = 5.66, 6.0) lattice for the 1 HYP overlap operator. Mp 
and fp are obtained from a fit to the (^o^o) correlator in the region 
T = {5Nt/16, . . . ,11Nt/1Q}- The solid curve is a linear fit, the dotted 
vertical lines indicate the physical M-^-ro and M^ro, respectively. 

(as shown in Fig.OJ as well as for m^ro and ms/rriud- Our final result is 

m,(MS, 2 GeV)ro = 0.301(25)(7) (11) 

in the continuum [x^/d.o.f. = 2.1] which, upon using tq = 0.5 fm, leads to 
the result given in @. For the ratio ms/rriud we find the value quoted in 
© [xVd.o.f. = 0.32]. In either case the generic comment on statistical and 
systematic errors, as stated above, applies. Evidently, any result involving 
mud benefits from our choice to use (|1()|) at finite lattice spacing, since it 
restricts the curve in Fig.|2to go through zero. The absence of curvature in 
our data is the reason why this ratio is in rather good agreement with the 
XPT result ms/rriud = 24.4(1.5) |2H|) in spite of the limitations mentioned 
and in spite of the latter value referring to a different theory (full QCD). 

The pseudoscalar decay constant fp is extracted from the (AqAq) cor- 
relator. Generically, we see plateaus for the effective decay constant from 
about Nt/4: on, again with the exception of the unfiltered operator where 
they are less pronounced. We use the Za values discussed earlier. Our 
central values stem from the interval 5A'7^/16 <t< IINt/IQ and the error 
estimate comes from comparing to the {PP) — {SS) channel, from varying 
the fit range and from the error on Z^. 

In Fig.|3] fp is plotted versus the meson mass squared. We see a strictly 
linear behavior at all couplings, with the lattice spacing having only mild 
effects on intercept and slope. Since (quenched) XPT predicts a linear de- 
pendence on Mp at lowest order and (quenched) chiral logs enter at 1-loop 
level only, we decided to stick to the leading order and do a linear fit. The 
intercept with {Mxrof = 1.251^ and {M^rof = 0.3537^ defines the fxro 
and /^ro at finite lattice spacing reported in Tabs.lHlandlTl respectively. 





(3 = 5.66 


(3 = 5.76 


(3 = 5.84 


(3 = 6.0 


none 
lAPE 
3 APE 
IHYP 
3 HYP 


ill-def. 
0.412(23)(69) 
0.426(18)(54) 
0.414(15)(40) 
0.400(15)(49) 


ill-def. 
0.391(21)(10) 
0.398(20)(06) 
0.398(21)(03) 
0.378(16) (08) 


0.575(39)(163) 
0.463(31)(35) 
0.449(33)(33) 
0.446(32)(22) 
0.419(28)(26) 


0.422(19)(10) 
0.430(22)(26) 
0.433(22)(31) 
0.430(22)(27) 
0.430(22)(35) 



Table 6: fK^o on all lattices used in the continuum extrapolation and the 
unfiltered. With the /?= 1 thin link action at /? = 5.66, 5.76 there is no signal. 





13 = 5.66 


(3 = 5.76 


13 = 5.84 


(3 = 6.0 


none 
1 APE 
3 APE 
IHYP 
3 HYP 


ill-def. 
0.367(37)(161) 
0.398(36)(133) 
0.373(18)(69) 
0.386(26)(118) 


ill-def. 
0.356(24)(22) 
0.358(24)(28) 
0.365(24)(14) 
0.334(18)(07) 


0.541(48)(211) 
0.438(45)(81) 
0.415(47)(66) 
0.412(42)(33) 
0.386(39)(55) 


0.387(23)(19) 
0.393(25) (36) 
0.397(26)(42) 
0.394(25)(36) 
0.396(25)(45) 



Table 7: f-j^VQ on all lattices used in the continuum extrapolation and the 
unfiltered. With the /9= 1 thin link action at /3 = 5.66, 5.76 there is no signal. 



0.45 




0.35 



Figure 5: Pseudoscalar decay constants f^^rQ and fRfo versus {a/rof. Errors 
are statistical, the 4 filterings are correlated (same configs.). 
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To extrapolate to the continuum, the same recipe is used as before. 
A combined fit to ah filtering levels with a common continuum limit and 
individual O(a^) terms is shown in Fig. |3 for fx and We obtain 

fKTo = 0.430(24)(5) (12) 

in the continuum [x^d.o.f. =0.96] which, via ro = 0.5fm, yields Q. Apply- 
ing the same strategy to fK/f-K we find the value quoted in Q [x^/d.o.f.= 
0.52]. Given the large errors it is not so surprising that both fx and fx/ fn 
are in reasonable agreement with the experimental f^"^ = 159.8(1.5) MeV 
and /x^//^^'' = 1.223(12), respectively [SHI. Especially in the latter case, 
higher precision would likely reveal a genuine quenching error [3Uj . 

4 Conclusion 

We have presented an exploratory scaling study with UV-filtered overlap 
quarks, concentrating on the quenched continuum limit of the light quark 
masses and the vr and K decay constants. The focus is not so much on the 
numerical values obtained; the point is that a nice scaling behavior has been 
observed down to lattices as coarse as ~ 1 GeV. As anticipated in the 
details of the filtering recipe play a minor role and therefore we conjecture 
that with analytic gauge link smearing techniques [HJ these properties will 
carry over to dynamical overlap calculations (cf . (HH 031 03 ) • 

For the time being let us mention that the simulations were done on a 
few of todays standard PCs (3 GHz Pentium 4), using an equivalent of ~120 
days on one such PC for the 1 APE operator (including all couplings) and 
of ~60 days for any of the other filtered operators. With a complete scaling 
study based on overlap valence quarks requiring less than 1 PC year one 
is tempted to ask whether for some phenomenological questions the overlap 
might actually represent the cheapest approach to obtain a continuum result. 

Acknowledgments: This work is based on an earlier collaboration with 
Urs Wenger that lead to [H]. S.D. was supported by the Swiss NSF. 
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